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1. Introduction 

Along the last two decades, the next-to-leading (NLO) corrections to the rapidity evolution 
of QCD scattering amplitudes in the asymptotic high-energy regime, given by the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) [1] equation and its nonlinear extensions, the Balitsky-Kovchegov (BK) 
[2] and JaliIian-Marian-Iancu-McLcrran-Wcigcrt-Lconidov-Kovner (JIMWLK) [3] equations, have 
become available [4-7]. Going beyond leading logarithmic accuracy is known to be mandatory in 
order to have a chance to accurately describe the available precision data on DIS at small Bjorken-x 
or forward pA collisions (for a review, see e.g. [8,9]). 

Flowever, it is known that the size of the NLO corrections to the BFKL kernel is very large, 
threatening to deprive the formalism of its predictiveness [10]. The source for this poor conver¬ 
gence was identified as due to large transverse col I incar logarithms appearing at each order in 
the perturbative expansion starting at NLO, and successful all-orders resummation procedures in 
Mellin space were devised [11] by effectively adding the higher-order terms involving poles in the 
characteristic function required by DGLAP evolution in collinear and anticollinear regions 1 . The 
convergence problem is not expected to disappear when including nonlinear effects (multiple scat¬ 
tering and saturation) as encoded in the BK equation [13], a fact that was confirmed explicitly in a 
recent numerical study [14,15]. 

Implementing the collinear improvement at the level of the non-linear BK equation is the major 
step earned out in the work [15], For this it is needed to identify the origin of the large corrections 
directly in coordinate space. Previous works in this direction [16,17] had already recognized the 
interplay between the kinematical constraint [18] and double collinear logs. In [15], we derived this 
relationship by a thorough analysis of Feynman diagrams in the double logarithmic approximation 
(Sec. 2). Although time ordering naturally provides a non-local evolution equation (2.8), we were 
able to show that rather nontrivially the double-log resummation can be implemented at the level 
of the (energy-independent) kernel, rendering its applicability more feasible than that of former 
approaches (Sec. 3). Flaving reformulated the resummation in terms of a modified BK kernel, we 
were able to show numerically its big impact on both stabilizing and slowing down the evolution 
(Sec. 4). This entails important consequences for phenomenology: in a recent analysis we show 
how resummed BK evolution provides an excellent fit to the combined HERA data for the DIS 
structure functions at small XBjorken with a small number of fitting parameters having values in the 
physically expected ballpark [19]. 

2. Double Logarithmic Approximation in High-Energy QCD Evolution 

The probability for gluon emission with transverse momentum k± and a fraction x of the parent 

parton longitudinal momentum is enhanced for soft and collinear emissions: && ~ cc^-rr- = 

k l 

cqdFdp, where a s = . Y = ln(l/x) is the rapidity variable and p = ln(/q /Ql), with Q o the 

’in the work [12], this procedure was approximately reformulated as an improvement at the level of the BFKL 
kernel in momentum space. 
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Figure 1 : Sample diagrams with two gluon ordered in longitudinal momentum (p + k 1 ) and also in 
lifetime (z p > Tj.): real-real graph (a); virtual-real graph (b). 


characteristic transverse scale of the target. In the kinematic region where a s Yp > 1 while a s Y <C 1 
and a s p <C 1, the resummation of terms of the form (a s Yp)" to all orders is expected to provide 
a reasonable approximation. This naive double-logarithmic approximation (DLA) appears in the 
BFKL equation, describing the rapidity evolution of the dipole S-matrix S x y = ^iv[V^ V y \ = \-T X y 

dyT xy (Y) = — J d Z-#ry Z [7i z (T) +T zy (Y ) — Txy(Y)\\ ^YY xyz = (2-1) 

by taking the limit where the daughter dipoles are much larger than the original one, \x — z\ — 
\z —y r = \x-y\. In this limit, averaging over angles and impact parameter T xy (Y) —> T(Y,r 2 ) 
and making the rescaling T(Y,r 2 ) = r 2 Q^s/ ( Y,r 2 ), Eq. (2.1) reduces to 

r Y r l /Qo d^ 2 

£/(Y,r 2 ) = £Y(0,r 2 ) + a s d yJ -Y-fJ(Y u z 2 ), (2.2) 

J0 Jr 2 

which resums by iteration the terms enhanced by factors ( a s Yp) n . A crucial point to notice, how¬ 
ever, is that these are not the only terms enhanced with the maximum number of large logs per 
power of the coupling. In fact, we will see that kinematical constraints also generate terms like 
a s p 2 ,a 2 Yp 2 ,a 2 p 4 , ■ ■ ■ which can be seen as large corrections to the impact factor or formally 
higher-order corrections to the BFKL kernel, but will become important and need to be resummed 
when Y > p > 1 / y/W s . To see how these higher powers of the transverse logs are generated, we 
consider a diagram like (a) in Fig. 1. It is proportional to 




d k+ 

~k+, 


dp 4 


'k+ P + p+ -)- k + £ p+ + k+ 


(p-k ) 2 

i 2 


Y =ln{q + /q^). 


(2-3) 


The factor +P 2 ' s ec l ua l to % + Tt ’ where z p = l/p = 2p + / p 2 is the associated gluon lifetime. 

To double log accuracy this can be replaced simply by 0( r p — T/ f ), and moreover this time ordering 
can be enforced directly in coordinate space [15]. Then, the DLA approximation to the sum of 
real-real graphs for NLO color dipole evolution reads 


CL 
2 K 




d k + 

~k+. 


dp 4 


'k+ 


’uz 


f'xyu 


uyz^xu^uz^zy 


X uzS xz S Z uS uy \e(p + a 2 - k + z ), (2.4) 


with u = max(|n — jc|, \u — y|); z = max(|z — jc|, \z — y|). As we saw in arriving to (2.2), the region 


z — x\ ~ \z — y| — |z — u\ ^$> \u — x\ ~ u — y| S> r = |jc — y| 


(2.5) 
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Figure 2: Anti-time-ordered diagrams ( T p < T&). 


generates the double logs. Indeed, within this approximation yz ^xyu 
T uz + T zy ~ 2 T(z 2 * ), which gives rise to the double logs through 


and 1 —S xu S uz S zy ~ 


r z " 2 dn 2 


r ci+ dp H 

A+4 p + 


d/r 


In^-ln^ 

u 1 


k+ 


= Fp-^; F = ln|pp = In^. (2.6) 


A quite non-trivial fact is that a number of anti-time ordered graphs (Fig. 2), involving factors 
p - +k ~ — 0(t* — T p ), which would also be potentially enhanced by double transverse logs since 


« 2 dd 2 

n <7 2 


dp 


+ 


'k+t P H 


fz 2 dn 2 ^ z 2 _ p 2 


2 ’ 


(2.7) 


cancel completely among themselves 2 . Thus we conclude that perturbative corrections enhanced 
by double logarithms Yp or p 2 can be resummed to all orders by solving a modified DLA equation 
involving manifest time-ordering 


£f(q + ,r 2 ) 


£/(0 ,r 2 ) + a s 



ho 


'? dlc + . , 

1F^<*' z 


( 2 . 8 ) 


3. The Resummed BK Evolution Kernel 


Eq. (2.8) is manifestly non-local in rapidity, being equivalent to dysY (F, p) = a s Jq dpi .c/(F — 
p +pi,p). Nevertheless, direct iteration of this equation and some mathematical manipulations 
allows us to write [15] 

*?(Y,p) = [ P dpJ(Y,p- Pl )*/(0,pi), f(Y,p) = f iZ r /( 7 ) e [«^A(r)i'+( 1 -/)P]. (3.1) 

Jo 2TU 

Here sY(Y,p) is the analytical continuation of the physical amplitude sY (F. p) for F > p, and 


OWfDLA (1") 


—(! — y) + \/ (! — 7) 2 + 4 ^ 


/(/) = l-a^DLA(7)- 


(3.2) 


2 Although virtual contributions do not involve a natural time ordering and are not enhanced by double logs, it is 

convenient to artificially split them in time- and anti-time-ordered contributions. In this way they regulate UV and IR 

divergences for both time orderings independently. 
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Figure 3: Numerical solutions to BK equation for the dipole amplitude {a s = 0.25,///((). p) = 1) at strict LO 
(left), including the NLO double log (center), and the resummed tower of double logs to all orders (right). 




Figure 4: Rapidity dependence of the target saturation momentum Q 2 S (Y), with a s = 0.25. 


Eq. (3.1) implies that we can write the DLA-resummed evolution of the amplitude in terms of a Y- 
independent kernel P£dla(p) (the inverse Mellin transform of the characteristic function 

A(Y,p) =A(0,p) + a s f d7i f dpi^D L A(p-pi)A(Fi,pi), ^dla(p) = ^\ (3.3) 
Jo Jo \/UsP 2 

Notice that the resummation also enters the initial condition, given by the limit of (3.1) at the an¬ 
alytically continued point Y = 0, /(0, p) = 8 (p) — \/W s J\ (2 \J a s p 2 ), from which analytical expres¬ 
sions for A (0, p) can be obtained for initial conditions of the Golec-Biernat-Wusthoff (&4 (0, p )gbw ~ 
1) [20] and McLerran-Venugopalan (j2^(0,p)mv ~ P) [21] type- After having performed the re¬ 
summation at double-logarithmic accuracy, it is easy to promote Eq. (3.3) into a more complete 
equation matching NLO BK 3 [15]: 

dyTxy = — j 

4. Some Numerical Results 

The numerical solution to Eq. (3.4) is presented in Fig. 3. We clearly see that while introducing 
the double log term from the NLO BK kernel renders the evolution unstable and physically mean¬ 
ingless, the evolution is perfectly smooth after resummation and certainly slower than that observed 
at LO. To characterize more precisely this reduction, we computed numerically the target saturation 

3 Single-log (collinear and running-coupling) and finite terms appearing in NLO BK have not been included in Eq. 
(3.4), although they can be simply added to the kernel. Prescriptions to resum single logs are discussed in [19]. 


(■ *-z ) 2 ] n (y-z) 


In 


(x-yp (x~y) 


[Txz + T Z y T xy T xz T zy ] 


Y > p. (3.4) 
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momentum Q%(Y) (Fig. 4). The growth of the saturation scale with Y is considerably reduced by 
the resummation: for sufficiently large rapidities, the saturation exponent X s = &{\n Q 2 S [Y)/Q^)/&Y 
is roughly reduced by a factor of 2 for the resunmied kernel as compared to LO one. This has 
important phenomenological consequences that have begun to be explored in [19], where also the 
relevant issue of single logarithmic (running coupling and collinear) corrections is addressed. 

References 

[1] L.N. Lipatov, Sov.J.Nucl.Phys. 23 (1976) 338; E.A. Kuraev, L.N. Lipatov, and V.S. Fadin, Sov.Phys. 
JETP 45 (1977) 199; I.I. Balitsky and L.N. Lipatov, Sov.J.Nucl.Phys. 28 (1978) 822. 

[2] I.I. Balitsky, Nucl.Phys. B463 (1996) 99; Yu.V. Kovchegov, Phys.Rev. D60 (1999) 034008. 

[3] J. Jalilian-Marian, A. Kovner, A. Leonidov, and H. Weigert, Nucl.Phys. B504 (1997) 415; Phys.Rev. 
D59 (1998) 014014; A. Kovner, J.G. Milhano, and H. Weigert, Phys.Rev. D62 (2000) 114005; E. 
Iancu, A. Leonidov, and L.D. McLennan, Nucl.Phys. A692 (2001) 583; E. Iancu, A. Leonidov, and 
L.D. McLerran, Phys.Lett. B510 (2001) 133; E. Ferreiro, E. Iancu, A. Leonidov, and L.D. McLerran, 
Nucl.Phys. A703 (2002) 489. 

[4] V.S. Fadin and L.N. Lipatov, Phys.Lett. B429 (1998) 127; M. Ciafaloni and G. Camici, Phys.Lett. 
B430 (1998) 349; V.S. Fadin (1998), [hep-ph/9807528]. 

[5] I.I. Balitsky, Phys.Rev. D75 (2007) 014001. 

[6] I.I. Balitsky and G.A. Chirilli, Phys.Rev. D77 (2008) 014019. 

[7] I.I. Balitsky and G.A. Chirilli, Phys.Rev. D88 (2013) 111501; A. Kovner, M. Lublinsky, and Y. 
Mulian, Phys.Rev. D89 (2014) 061704. 

[8] Yu.V. Kovchegov and E. Levin, Quantum Chromodynamics at High Energy, Cambridge University 
Press, Cambridge 2012. 

[9] F. Gelis, E. Iancu, J. Jalilian-Marian, and R. Venugopalan, Ann.Rev.Nucl.Part.Sci. 60 (2010) 463. 

[10] D.A. Ross, Phys.Lett. B431 (1998) 161. 

[11] G. Salam, JHEP 9807 (1998) 019; M. Ciafaloni, D. Colferai, and G. Salam, Phys.Rev. D60 (1999) 
114036; M. Ciafaloni, D. Colferai, G. Salam, and A.M. Stasto, Phys.Rev. D68 (2003) 114003. 

[12] A. Sabio Vera, Nucl.Phys. B722 (2005) 65. 

[13] E. Avsar, A.M. Stasto, D.N. Triantafyllopoulos, and D. Zaslavsky, JHEP 1110 (2011) 138. 

[14] T. Lappi and H. Mantysaari, Phys.Rev. D91 (2015) 074016. 

[15] E. Iancu, J.D. Madrigal, A.H. Mueller, G. Soyez, and D.N. Triantafyllopoulos, Phys.Lett. B744 
(2015) 93. 

[16] L. Motyka and A.M. Stasto, Phys.Rev. D79 (2009) 085016. 

[17] G. Beuf, Phys.Rev. D89 (2014) 074039. 

[18] M. Ciafaloni, Nucl.Phys. B296 (1988) 49; B. Andersson, G. Gustafson, and J. Samuelsson, 

Nucl.Phys. B467 (1996) 443; J. Kwieciriski, A.D. Martin, and P.J. Sutton, Z.Phys. C71 (1996) 585. 

[19] E. Iancu, J.D. Madrigal, A.H. Mueller, G. Soyez, and D. N. Triantafyllopoulos, Collinearly-Improved 
BK Evolution Meets the HERA Data [arXiv: 1507.03651 [hep-ph] ]. 

[20] K.J. Golec-Biernat and M. Wusthoff, Phys.Rev. D59 (1998) 014017. 

[21] L.D. McLerran and R. Venugopalan, Phys.Rev. D49 (1994) 3352. 


6 



